Rotational tuning of interaction in metamaterials 
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We experimentally observe the tuning of metamaterials through the relative rotation of the elements about 
their common axis. In contrast to previous results we observe a crossing of resonances, where the symmetric 
and anti-symmetric modes become degenerate. We associate this effect with an interplay between the magnetic 
and electric near-field interactions and verify this by calculations based on the interaction energy between 
resonators. 
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Metamaterials created as an array of sub-wavelength, 
resonant elements can exhibit interesting electromagnetic 
properties, such as a negative refractive index-. Unlike 
the atoms in natural materials, the near-field patterns 
of metamaterial elements are quite complex, giving rise 
to strong interactions between them. Understanding this 
interaction is essential as it determines the overall reso- 
nant properties and effective parameters of the material. 
By controlling the relative arrangement of elements, it is 
possible to change this coupling and tune the properties 
of the structure's. This allows us to alter the response 
of the material substantially without having to greatly 
alter the geometry or constituents. 

An important building block of metamaterials is the 
Split Ring Resonator (SRR) 5 , which can exhibit a nega- 
tive magnetic response due to its strongly resonant mag- 
netic polarizability. By coupling a pair of SRRs, chiral 
properties^ and trapped dark mode*" can also be ob- 
served. In this Letter, we study the dynamics of two 
microwave SRRs broadside coupled to each other, vary- 
ing twist angle 9 between them, as shown in Fig. [lja). 
A similar system operating at near infra-red frequencies 
was presented in Ref. [21 however in this work we show 
a substantially different regime of interaction verified ex- 
perimentally for microwaves. 
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FIG. 1. (a) A schematic showing the rings rotated with re- 
spect to each other through angle 0, and the polarization of 
the incoming waves, (b) A comparison of the experimental 
(solid) and numerical (dashed) absorption for angle 6 = 90°. 



The rings used have an inner radius of 3.5mm, an outer 
radius of 4mm, and a gap of 1mm. They are copper, 
printed onto 1.6mm thick FR4 circuit board, and the 
rings are 3.6mm apart, with the dielectric boards located 
between the rings. The incoming microwaves are polar- 
ized so that the electric field is in the x direction, as 
shown in Fig. [l|a). Numerical calculations were per- 
formed using CST Microwave Studio with the boards 
having a dielectric constant of 4.6. 

The most meaningful definition of the resonant fre- 
quency is the frequency of maximum excitation of cur- 
rents within the rings. This can be determined most 
readily by considering the absorption of the system, given 
by 1— l^il 2 — |£ii| 2 , where S21 is the transmission coeffi- 
cient, and S11 is the reflection coefficient. A comparison 
of the experimental and numerical curves for 9 = 90° is 
shown in Fig. WCb) . 

Experimental results were measured using a Rohde and 
Schwarz ZVB network analyzer in a WR-229 rectangu- 
lar waveguide, with 9 varied from 0° to 180° in 10° in- 
crements, while numerical results were calculated in 5° 
increments. The resonant frequency for each angle was 
found from the maximum of the absorption curve, and 
the resulting numerical and experimental results are com- 
pared in Fig. [2j 
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a ' Electronic mail: david.a.powell@anu.edu.au 
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FIG. 2. A comparison of the experimental (solid line) and 
numerical (dashed line) resonant frequencies. 
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FIG. 3. (a) The numerical absorption curve for 6 C (33°), and 
(b) the numerical (dashed) and experimental (solid) absorp- 
tion curves for 9 = 180°. 



For 9 — 0° , there are two resonances us and ujaSi an d 
by inspection of the currents in the rings we verify that 
these correspond to the expected symmetric and anti- 
symmetric modes. As 9 increases, loas increases and uis 
decreases, reaching their maximum and minimum values 
respectively at 6 = 180°. For our chosen parameters 
the resonances cross at 9 C s» 33°, in contrast to Ref. 2, 
where an avoided crossing of resonances was found which 
was attributed to the electric quadrupole and octupole 
moments of the rings. 

The numerical absorption curve at 33° is shown in 
Fig.ga). This curve was calculated assuming loss free 
dielectric boards, so as to ensure the position of the cross- 
ing, given the low coupling of the asymmetric mode close 
to 9 C . It is quite clear here that there is only one reso- 
nance to be observed. 

As can be seen in Fig. [31(b) , the two resonances have 
different Q factors (widths), which can be attributed to 
conduction and radiation losses. The suppression of ra- 
diation losses in the symmetric mode is related to the 
magnetic and electric dipole arrangements^, where the 
electric dipole has much stronger radiation losses. The 
asymmetric mode has mostly an electric dipole response, 
causing greater radiation losses, whereas the radiation 
losses at the symmetric mode have been suppressed by a 
dominating magnetic dipole. The radiation distribution 
then changes with angle 9, as can be seen by comparing 
Figs. Iljb) and §b). 

The tuning of the system by rotation can be explained 
by looking at the interaction between the rings. As the 
rings are twisted, the magnetic and electric near-fields of 
the two rings change, changing the coupling between the 
resonances. This is approached theoretically using the 
Lagrangian for a pair of coupled SRRS 

C = A{Q\ + Ql + 2aQ 1 Q 2 )~B{Ql + Ql + 2!3Q l Q 2 ) (1) 

where a and j3 are the dimensionless magnetic and elec- 



tric interaction constants and Q(t) is the time-dependent 
mode amplitude. By substituting Eq. ([I]) into the Euler- 
Lagrange equation we find that 



lo Qi = 

, ,2/ 



-aQ 2 



(2) 
(3) 



Qi 

Q 2 + ujqQz = -otQx - P^oQi 

This then allows the two resonances to be found - sym- 
metric (when Qi = Q 2 ), and asymmetric (when Q\ = 
-Q 2 ): 
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In principle if ojq is known, then by inverting Eq. Q it 
is possible to fit a and /3 from u>s and u>as- However we 
found that this procedure is extremely sensitive to error 
and does not yield usable results. Instead, we start from 
the approach outlined in Ref. 4]to evaluate the interaction 
energy between the fundamental modes of the rings. 

For a pair of rings in a homogeneous dielectric back- 
ground, the interaction constants are shown in Fig. 4][a). 
We find that the interaction constants are very well de- 
scribed by P = Pi cos(#) and a — a + ct\ cos((9) with 
/3i = 0.085, a = 0.098 and a x = 0.05. These constants 
are dictated by the charge separation across the gap of 
the ring, the current circulating around the ring, and 
the inhomogeneity of the current distribution around the 
ring, respectively. For rings aligned on the same axis, we 
expect that the magnetic interaction should always be 
positive, as the intersecting magnetic field from one loop 
should always be normal to the other loop. In addition 
the electric interaction should be positive at 9 = 0° as 
the charge distribution has the nature of parallel dipoles. 
All arrangement of rings on the same axis which we con- 
sidered obeyed these considerations. 

In Fig. Qb) we plot the corresponding frequencies of 
the symmetric and anti-symmetric modes, normalized to 
loq. As our approach models the response of the res- 
onators in a homogeneous dielectric background, the re- 
sults are significantly different from those observed exper- 
imentally, where the dielectric is inhomogeneous and the 
effect of waveguide boundaries is also significant. In par- 
ticular, the crossing of resonances cannot be reproduced 
for rings in homogeneous background. Therefore we con- 
sider the possible regimes of interaction which may occur, 
under the assumption that the interaction constants can 
be fitted as described above. 

The case considered in Fig. (4^a-b) corresponds to the 
magnetic interaction always being larger than the elec- 
tric interaction. This results in increasing splitting of ujg 
and luas with increasing twist angle, however in prin- 
ciple there is no reason why the splitting cannot de- 
crease. We show such a case in Fig. 4|c-d), where we 
have set /?i = 0.02 < ai, such that the inhomogeneity 
in the current has a stronger influence than the dipole- 
like charge distribution. Despite the apparent difference 
in frequency splitting curves, there is little difference be- 
tween the interaction constants shown in Fig. Ba) and 
(c). 
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FIG. 4. (a) Magnetic (a) and electric (/3) interaction con- 
stants calculated for a pair of rings in free space, and (b) 
corresponding resonant frequencies, (c) Interaction constants 
for case where electric coupling dominates at 9 — 0°, and (d) 
corresponding resonant frequencies, (e) Interaction constants 
which become equal for some angle 8, and (f) corresponding 
resonant frequencies. 



not justified. However by taking the inhomogeneity of 
magnetic interaction into account, we find that we can 
neglect higher order electric interactions. 

Equations Q show that the tuning curves arise from 
competition between electric and magnetic interaction 
constants. This means that even when interaction is 
strong, the corresponding frequency splitting can still 
be weak. The experimentally-observed crossing angle 
8 C = 33° represents a particularly interesting case, and 
it occurs when a = (3. In the lossless case the symmetric 
and anti-symmetric resonances become degenerate, and 
although both electric and magnetic interaction coeffi- 
cients are strong, they effectively cancel each other out. 
However in the presence of losses, we need to consider 
the eigenfrequencies lo$ and loas as complex values, with 
imaginary parts corresponding to the losses. By exam- 
ining Fig. |3jb) we see that the coupling of the modes to 
the waveguide is not identical, hence they have different 
radiation losses. Thus it becomes very unlikely that in 
the complex plane lo$ — loas for any parameter value, 
and true degeneracy does not occur. Therefore it is not 
possible to cancel out the interaction between rings, and 
our numerical work verified this. 

In conclusion, we have shown that by changing the 
relative rotation between two rings, we can significantly 
change the coupling, which causes the resonances to 
change. We have found that there is a crossing where 
the two resonances coexist, which corresponds to equal 
electric and magnetic coupling. 
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The only other case allowed in our model of interaction 
under the afore-mentioned physical constraints on a and 
f3 is that a > j3 for 9 = 0° . An example of this is given in 
Fig. Qe), where we have reduced the magnetic coupling 
such that at some angle a = (3, by setting ao = 0.04, a% = 
0.02. The corresponding resonant frequencies normalized 
to loq are plotted in Fig. Qf ) . We see that for low twist 
angle, lo$ occurs at a higher frequency, but decreases with 
angle and crosses loas- Clearly this regime corresponds 
to what we observe in experiment, and we hypothesize 
that the inhomogeneous dielectric serves to enhance the 
electric interaction relative to the magnetic interaction 
in our system. 

In Ref. [2] it was claimed that the resonances converge, 
undergo an avoided crossing, then diverge, as 9 increases. 
Utilizing both our analytical model and full numerical 
simulation, we could not find this regime in our sys- 
tem. We also note that for our structure the assump- 
tion of constant magnetic interaction with twist angle is 
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